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ABSTRACT 


The  flow  field  of  nozzles  exhausting  into  a  vacuum  was  investigated 
both  analytically  and  experimentally.  The  applicability  of  the  method 
of  characteristics  for  jet  expansions  into  a  free  molecular  environment 
is  discussed.  Fifteen- degree,  half-angle,  conical  nozzles  with  area 
ratios  of  1.  0  to  207  were  tested  with  argon  and  nitrogen  at  total  pres¬ 
sures  of  100,  50,  and  35  psia.  Ambient  pressures  of  10"^  to  10"®  torr 
yielded  total  pressure  to  ambient  pressure  ratios  of  10^  to  10^.  Prop¬ 
erties  of  the  flow  field  were  determined  by  the  use  of  flat  plates  and  flow 
direction  vanes. 
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NOMENCLATURE 

A  Area 

A+  Nozzle  throat  area 

Cd  Drag  coefficient 

F  Normal  force  on  the  flat  plate 

K.n  Knudsen  number  based  on  flat  plate  dimension 

K. n “  Knudsen  number  based  on  throat  radius 

M  Mach  number 

M  Mass  flow  rate 

F  Pressure 

It  Gas  constant 

i'  Radial  distance  from  axis 

r*  Throat  radius 

T  Temperature 

V  Velocitv 

X  Axial  distance  from  nozzle  exit 

y  Ratio  of  specific  heats 

,5  Boundary  layer  thickness 

Boundary  layer  displacement  thickness 
0N  Nozzle  half  angle 

\  Mean  free  path 

v  Prandtl- Meyer  angle 

p  Density 

Nozzle  chamber  conditions 
Nozzle  exit  conditions 
Maximum  conditions 

Standard  temperature  pressure  and  density  {300°K, 
1  amagat) 

Background  conditions 


SUBSCRIPTS 


c 

Max 


Vlll 


1  atm. 
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SECTION  I 
INTRODUCTION 


When  the  exhaust  from  a  nozzle  discharges  into  a  low  ambient  pres¬ 
sure,  it  expands  to  a  much  larger  extent  than  at  sea  level  pressure.  It 
is  desirable  to  study  this  plume  for  several  reasons.  The  effects  of  the 
jet  plume  on  a  vehicle  system,  such  as  impingement  on  vehicle  sur¬ 
faces  and  radio  wave  attenuation,  require  knowledge  of  the  plume  prop¬ 
erties.  The  use  of  exhaust  scavenging  systems  in  space  environmental 
chambers  requires  detailed  information  of  the  jet  properties.  Also  the 
free  expanding  jet  may  be  used  as  a  research  tool  for  studying  the  rare¬ 
fied  gas  flow  regime  if  its  properties  can  be  determined. 

Most  of  the  past  measurements  (Ref.  1,  2,  and  others)  have  been 
made  to  determine  the  shape  of  the  boundary  of  the  jet  by  observing  the 
boundary  shock  wave  by  optical  methods.  However,  at  the  pressures 
characteristic  of  the  space  environment,  the  mean  free  path  is  much 
too  great  for  a  shock  wave  to  be  defined.  Thus,  it  becomes  necessary 
to  survey  the  jet  flow  field  to  determine  the  spatial  distribution  of  its 
properties,  i.  e. ,  velocity,  density,  etc. 

The  present  study  is  an  extension  of  a  development  test  for  the 
Nimbus  attitude  control  jet  reported  in  Ref.  3.  In  that  work,  the  jet 
from  a  specific  nozzle  was  surveyed  by  determining  the  normal  force 
on  a  small  flat  plate.  The  tests  reported  herein  are  primarily  con¬ 
cerned  with  determining  the  applicability  of  the  method  of  characteris¬ 
tics  for  predicting  the  flow  field  and  developing  techniques  for  rarefied 
flow  investigations. 


SECTION  II 

CALCULATION  OF  THE  FLOW  FIELD  DOWNSTREAM  FROM  A 
NOZZLE  EXHAUSTING  INTO  A  LOW  PRESSURE  ENVIRONMENT 


2.1  METHOD  OF  CHARACTERISTICS 

The  method  of  characteristics  is  a  computational  procedure  by 
which  the  properties  of  a  flow  field  of  a  compressible  fluid  can  be 
determined.  The  state  of  the  ideal  exhaust  gas  and  its  local  velocity 
vector  can  be  mapped  over  a  field  defined  by  the  boundary  conditions. 
In  essence,  a  point  in  the  field  has  a  zone  of  influence  over  a  region 
downstream  within  a  wave  which  makes  an  angle  sin-1  (acoustic  speed/ 
stream  speed)  with  the  velocity  vector.  Through  the  isentropic  flow 
relationships,  the  effects  on  flow  velocity  and  gas  state  within  the  zone 
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of  influence  can  be  evaluated  (Ref.  4  and  5).  The  boundaries  of  the 
zones  of  influence  are  termed  characteristic  lines.  The  behavior  of 
the  flow  properties  along  these  lines  is  governed  by  differential  equa¬ 
tions  which  are  put  in  finite  difference  form  for  computation.  The 
program  used  in  the  present  calculations  was  developed  for  the  IBM  7090 
computer  by  Lockheed  Missile  and  Space  Company  (Ref.  6).  The  com¬ 
putation  starts  with  uniform  flow  properties  on  a  spherical  surface  at  the 
nozzle  exit.  The  other  initial  condition  is  defined  by  a  Prandtl- Meyer 
expansion  of  the  flow  at  the  nozzle  lip  to  the  ambient  pressure,  Pw. 


-  Nozzle  Half-Angle 

v00-  Prandtl-Meyer  Angle  Corresponding  to  Background 
to  Stagnation  Pressure  Ratio 

ue  -  Prandll-Meyer  Angle  Corresponding  to  Nozzle  Exit 

Mach  Number  / 

r*  =  Throat  Radius 
Ue  =  Velocity  at  Exit 


From  a  number  of  points  at  the  exit,  a  network  of  characteristics 
lines  is  built  up  step-by-step  to  yield  the  flow  properties  at  the  inter¬ 
section  points.  The  results  of  the  computation  are  coordinates  of  the 
intersection  points,  Mach  number,  and  flow  inclination  to  the  axis.  Since 
the  flow  is  isentropic,  the  Mach  number  uniquely  relates  the  density, 
pressure,  temperature,  and  velocity  to  the  state  in  the  reservoir.  A  sub¬ 
program  devised  by  the  Scientific  Computing  Branch  of  AEDC  interpolated 
the  data  to  provide  contours  of  constant  Mach  number  and  constant  mass 
flow  stream  tube  boundaries.  Figure  1  shows  a  typical  example  of  the 
form  of  the  computed  results.  The  spread  of  the  computed  points  about 
the  faired  contour  results  from  a  finite  mesh  size.  The  numerical  diffi¬ 
culties  in  the  vicinity  of  the  axis  are  mentioned  by  Ferri  (Ref.  5)  and  re¬ 
quire  judgment  in  the  fairing  of  the  contours. 
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The  variables  of  the  problem  are:  ratio  of  ambient  pressure  to 
nozzle  chamber  pressure  ( PK / Pc )  ,  gas  specific  heats  ratio  iy),  nozzle 
area  ratio,  and  nozzle  exit  half-angle  (0n)  ,  For  the  present  calcula¬ 
tion,  the  nozzle  half-angle  was  constant  at  15°,  and  the  area  ratio  varied 
from  3.  73  to  207  to  conform  to  the  experiments.  Two  specific  heat 
ratios,  y  =  1.4  and  1.  67,  which  apply  to  diatomic  and  monatomic  gases 
were  used.  No  significant  difference  between  the  computed  flow  prop¬ 
erties  within  the  m/m*  =  .  9  contour  for  the  two  pressure  ratios 
( PM  /'Pc  =  2  x  10“e  and  2  x  10“ 7 )  was  observed,  and  it  can  be  concluded 
that  the  flow  field  comprising  at  least  90  percent  of  the  mass  flux  would 
be  the  same  if  the  ambient  pressure  were  reduced  to  zero. 

The  Mach  number  loses  its  physical  significance  in  a  rapidly  expand¬ 
ing  flow.  The  constant  Mach  number  contours  are  replaced  in  Figs.  2 
through  11  with  constant  density  ratio  curves  that  can  be  interpreted  as 
mass  or  number  density.  These  figures  give  a  description  of  the  momen¬ 
tum  flux  field,  which  determines  the  forces  and  heat  transfer  to  surfaces 
in  the  flow.  It  is  important  to  observe  that  beyond  200-300  throat  radii 
from  the  exit  the  streamlines  become  practically  ^traight  and  the  density 
variation  is  inverse  with  the  square  of  distance  (-.source  flow).  The 
source-like  flow  is  established  in  about  100  throat  radii  for  y  =  1.  67  and 
by  400  for  y  =  1.4.  If  the  flow  behaves  ideally  to  a  density  ratio  of 
(>/pc  =  1  *  TO-1,  then  subsequent  deviations  (such  as  non- continuum  or 
internal  mode  freezing  effects)  will  not  influence  the  flow  field  shown. 

The  field  can  therefore  be  extended  to  infinity  for  the  zero  ambient  pres¬ 
sure  situation. 

In  the  present  experiments  (PM  ~  I0~*  to  I0‘*a  iorr)  ,  the  mean  free 
path  in  the  ambient  region  is  about  the  length  of  the  flow  field  shown  in 
the  figures  and  about  the  same  as  along  the  p;pc  <  10“ 7  contour.  A 
boundary  shock,  characteristic  of  higher  ambient  pressures,  could  not 
form  and  the  interference  between  the  nozzle  flow  and  background  gas  is 
simply  a  linear  addition  in  density  with  no  significant  momentum  exchange. 

Figures  12  and  13  show  the  effect  of  nozzle  area  ratio  on  the  expan¬ 
sion  of  the  flow  field  for  diatomic  and  monatomic  gases.  As  more  ex¬ 
pansion  occurs  within  the  nozzle  -  i.e.,  higher  area  ratios  -  the  closer 
the  flux  will  be  to  the  axis.  The  effect  of  y  reflects  the  distribution  of 
energy  in  the  internal  modes  of  the  gas  molecular  motion.  At  the  nozzle 
exit,  the  density  ratio  is  about  the  same  for  the  diatomic  and  monatomic 
gas.  However,  the  temperature  of  the  monatomic  gas  is  only  60  percent 
of  the  other;  thus  it  is  much  nearer  to  full  expansion  (static  tempera¬ 
ture - *-0)  at  the  nozzle  exit,  and  the  jet  lies  closer  to  the  axis.  A  limita¬ 

tion  of  the  numerical  solution  by  the  method  of  characteristics  is  the 
convergence  of  the  characteristics  lines  at  high  expansions.  The  angle 
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between  the  characteristics  and  the  streamlines  is  the  Mach  angle 

(s‘n  1  =  8in~'  -j- )v-  -At  large  Mach  number,  the  intersections  of  the 

characteristic  net  become  far  apart  and  poorly  defined.  The  numerical 
difficulty  combined  with  the  fact  that  the  ideal  representation  of  the  gas  is 
no  longer  physically  realistic,  dictates  terminating  the  solution  at  a  Mach 
number  of  40  in  the  program.  This  is  the  reason  that  the  monatomic  flow 
fields  are  reproduced  to  a  larger  scale  than  are  the  diatomic  (y  =  1.4) 
fields. 

In  summary  the  method  of  characteristics  is  based  on  the  following 
s  impl  if  ic  at  ions : 

(a)  The  gas  has  an  ideal  isentropic  behavior,  i.  e.  ,  it  can  be  re¬ 
presented  as  having  a  constant  ratio  of  specific  heats  beyond 
the  nozzle  exit. 

(b)  The  gas  is  inviscid  or  the  effects  of  viscosity  are  negligible. 

<c)  The  gas  is  a  continuous  medium. 

The  departure  of  the  real  physical  situation  from  these  conditions 
is  discussed  in  the  following  sections: 

2.1.1  Condensation 

Condensation  may  occur  within  the  nozzle  or  beyond  the  exit  plane. 
Within  the  nozzle  it  can  take  the  form  of  a  condensation  shock  resulting 
from  the  release  of  the  heat  of  condensation.  The  Mach  number  behind 
the  shock  will  be  reduced  to  a  subsonic  or  a  lower  supersonic  value 
(Ref.  7)  with  the  result  that  the  nozzle  exit  Mach  number  is  less,  and 
jet  expansion  outside  the  nozzle  will  be  larger  than  predicted  by  the 
isentropic  calculation.  The  experimental  results  for  argon  strongly 
suggest  that  condensation  caused  the  large  expansion. 

The  condensation  of  nitrogen  is  discussed  in  detail  in  Appendix  I. 

It  is  concluded  that  condensation  of  nitrogen  in  the  nozzle  will  not  occur 
in  the  present  geometry  when  the  area  ratio  is  less  than  70.  For  the 
two  larger  area  ratios  (Ae/A*  =  70  and  200),  the  data  on  which  the  esti¬ 
mates  were  made  leave  uncertainty  as  to  the  probability  of  condensation. 
There  are  insufficient  data  for  argon  to  make  a  similar  prediction. 

Beyond  the  nozzle  exit,  the  collision  frequency  decreases  rapidly 
and  condensation  is  less  probable  than  in  the  nozzle.  However,  the  total 
number  of  collisions  is  proportional  to  size  of  the  nozzle  and  supply 
stagnation  density.  If  condensation  occurs  in  the  intermediate  field  be¬ 
tween  the  exit  and  the  source -like  outer  part,  a  two  phase  separation  of 
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flux  would  be  expected  with  the  gaseous  component  expanding  further  than 
for  no  condensation.  The  solid  part  of  the  mass  flux  would  be  more 
concentrated  near  the  axis. 

If  condensation  appears  in  the  field  where  the  streamlines  are 
already  linear,  it  would  be  impossible  to  detect  it  by  momentum  flux 
measuring  techniques  such  as  are  presently  employed.  The  energy  in 
modes  other  than  directed  velocity  is  so  small  that  condensation  would 
be  a  minor  perturbation. 


/  2.1.2  Viscosity 


Viscous  effects  on  the  flow  field  will  have  an  influence  where  velocity 
gradients  or  large  rates  of  change  of  density  are  present.  These  con¬ 
ditions  exist  at  the  nozzle  wall  and  at  the  exit  lip  where  the  corner  of  the 
Prandtl- Meyer  expansion  starts.  The  most  obvious  effect  of  the  nozzle 
boundary  layer  is  to  reduce  the  effective  area  ratio  below  the  geometrical 
area  ratio.  The  subsonic  portion  of  the  boundary  layer  will  expand  into 
the  background  gas  to  help  establish  the  ambient  pressure.  The  flow  out¬ 
side  the  boundary  layer  will  turn  through  the  Prandtl- Meyer  angle  deter¬ 
mined  by  exit  Mach  number  and  ambient  pressure.  The  supersonic  region 
of  the  boundary  layer  will  turn  through  an  angle  determined  by  the  local 
Mach  number.  The  lip  of  the  nozzle  will  be  effectively  rounded  by  the 
boundary-layer  (see  figure  below). 


The  primary  effect  of  the  nozzle  boundary  layer  is  to  move  the  limit 
expansion  line  downstream  from  the  lip  leaving  a  gap  through  which  the 
boundary  layer  gas  expands  into  the  ambient  background.  The  change 


( o 


5 


A  E  DC-T  R-65-15 1 


from  the  inviscid  flow  field  to  that  with  a  boundary  layer  is  indicated 
in  the  figure  below. 


If,  for  example,  the  boundary  layer  is  10  percent  of  the  total  flow, 
the  90  percent  mass  flux  streamline  takes  the  slope  of  the  limit  expansion 
line,  and  there  is  a  slight  displacement  of  the  other  mass  flux  lines.  The 
constant  density  lines  are  shifted  toward  the  axis  for  m/m*  >  0.8  since 
the  total  flux  between  80  percent  and  90  percent  remains  the  same.  The 
effect  of  the  nozzle  boundary  layer  will  be  to  reduce  the  momentum  flux 
near  the  outer  edge  of  the  inviscid  jet.  The  effective  area  ratio  of  the 
nozzle  will  be  reduced  from  the  geometrical  area  ratio  by  the  boundary 
layer  displacement  thickness.  The  nozzle  boundary  layer  has  been  com- 
puted  by  the  Cohen-Reshotko  method  (Ref.  8)  for  the  range  of  pressures 
used  in  the  experiments.  Figures  14a  and  14b  show  the  results  in  terms 
of  the  effect  on  nozzle  Mach  number  and  the  boundary  layer  thickness. 

At  the  same  supply  pressure,  boundary  layer  will  be  thicker  for  argon 
than  for  nitrogen. 

The  estimated  effect  of  the  boundary  layer  on  exit  Mach  number  is 
small  (<  4%)  for  nitrogen  and  would  not  produce  an  appreciable  influence 
on  the  flow  field.  However,  when  the  viscous  layer  reaches  10  percent 
of  the  nozzle  radius,  20  percent  of  the  mass  flow  will  be  entrained.  There¬ 
fore,  the  flow  field  outside  the  streamtube  for  0.  8  mass  flux  will  be  modi¬ 
fied  by  the  effects  of  viscosity.  This  is  the  situation  when  nitrogen  is  the 
gas  at  a  pressure  of  2. 4  atm  in  the  area  ratio  70  nozzle  and  about  5  atm 
for  the  largest  area  ratio  (207). 

2.1.3  Continuum  Limitations 

One  of  the  interesting  aspects  of  the  expansion  of  a  jet  into  a  vacuum 
is  that  the  gas  traverses  smoothly  through  a  series  of  gas  dynamic 
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regimes  that  are  usually  represented  by  different  analytical  models: 

(a)  the  continuum,  (b)  a  transition,  and  (c)  free-molecule.  These  are 
distinguished  by  the  magnitude  of  the  Knudsen  number:  the  ratio  of  the 
mean  free  path  to  some  significant  body  or  boundary  surface  dimen¬ 
sion.  If  the  distance  between  molecular  collisions  is  small  relative  to 
the  size  of  the  surface,  a  continuum  model--  such  as  the  method  oi 
characteristics- -is  applicable.  In  free-molecule  flow,  the  inter- 
molecular  collisions  have  no  influence  on  the  flow  about  the  body  in 
the  stream.  The  streamlines  are  linear  while  the  paths  of  the  individual 
molecules  are  the  resultant  of  the  mean  stream  velocity  and  the  random 
velocity  corresponding  to  the  translational  temperature. 

In  the  transition  regime,  there  is  an  interaction  of  one  region  of  the 
flow  field  on  another.  This  region  has  been  attacked  from  the  free- 
molecule  limit  for  deriving  pressures  on  bodies  immersed  in  the  flow. 

The  free-molecule  flow  is  modified  in  the  'first  collision"  theory  by 
impacts  of  surface  emitted  molecules  with  those  approaching  the  surface. 
However,  in  the  case  of  a  free  expansion  field,  it  is  difficult  to  designate 
a  significant  dimension  to  define  a  Knudsen  number.  The  radius  of  cur¬ 
vature  of  a  streamline  might  be  considered,  but  the  functional  relation¬ 
ship  between  such  an  arbitrary  dimension  and  the  mean  free  path  is  not 
obvious. 

There  is  a  close  analogy  between  the  disturbance  wave  of  the  con¬ 
tinuum  theory  and  the  particle  velocity  in  the  actual  gas  that  suggests 
that  the  transition  regime  is  not  sharply  distinguishable  for  an  expan¬ 
sion.  For  the  continuum,  the  zone  of  influence  is  a  cone  whose  half¬ 
angle  is  sin*1  for  the  kinetic  model  gas,  the  mean  zone  of  influence 

is  also  a  cone  about  the  velocity  vector  sin-1  .  “a”  is  the  acoustic 

speed,  \j y  and  u  is  the  mean  molecular  speed,  a  / -fi-  -Jbi-  ,  while  u 

is  the  flow  speed.  The  similarity  of  these  quantities  and  the  fact  that 
macroscopic  properties  of  the  gas  can  be  predicted  from  the  kinetic  molec¬ 
ular  model  lead  to  the  expectation  that  the  transition  from  a  continuum  flow 
model  at  the  nozzle  exit  to  a  free  molecule  state  downstream  will  be  nearly 
the  same  as  the  characteristics  method  results. 

A  conservative  criterion  for  applicability  of  the  continuum  analysis 
can  be  specified  as  follows:  Kn  =  .1  ,  for  a  100  r*-long  segment,  a  , 

of  the  contour,  =  I0“s.  Then  Kn  =  —  =  ■  10  2  Kc  =  1  03  Kn*  or 

KnK  ^Kn*  =  ^  is  10~4  .  For  the  test  nozzle  with  a  0.  5  mm  throat,  the 

reservoir  density  must  be  unity  or  greater  (the  product  r*  >  0.5 
amagat  -  mm. 
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2.1.4  Freezing  of  Internal  Energy  Modes 

The  internal  energy  of  a  gas  is  distributed  among  several  modes: 
ionization,  chemical  (commonly  dissociation  for  poly-atomic  mole¬ 
cules),  vibration,  rotation,  and  translation  of  the  molecules .  In  equilib¬ 
rium  these  modes  become  inactive  at  progressively  lower  temperatures 
and  influence  the  form  of  the  state  equation.  When  a  gas  undergoes 
a  rapid  expansion,  the  time  to  deactivate  the  various  modes  is  important, 
and  it  is  reasonable  to  associate  this  time  with  the  number  of  particle 
collisions.  There  is  approximately  one  collision  for  translation  and  10 
for  rotation,  etc.  It  is  usual  to  consider  the  expanding  gas  an  ideal 
(const,  y )  medium  until  "freezing"  of  a  mode  occurs.  At  this  point  a 
new  y  describes  the  expansion  process  until  the  next  mode  freezes. 


Vibrational  equilibrium  in  nozzles  is  considered  in  Ref.  9.  The 
nozzle  supply  temperature  was  near  300°K  for  the  present  experiments 
so  that  only  rotation  and  translation  energies  are  important  for  the  dia¬ 
tomic  gas  (N2)  ■  As  the  expansion  progresses,  the  collision  rate  drops 
(cf  Appendix  I),  and  equilibration  of  the  rotational  mode  ceases  to  be 
important.  The  gas  then  acts  like  a  monatomic  medium  with  y  =  1.  67 
rather  than  1.4.  Further  expansion  results  in  translational  freezing 
which  is  equivalent  to  free-molecule  flow.  Reference  10  derives  a 
dimensionless  parameter  to  predict  the  rotational  freezing  Mach  num¬ 
ber  in  an  isentropic  expansion  from  a  sonic  orifice: 


run 


JM 

il  (  x  /  d  *  J 


d*  is  the  throat  diameter,  a^is  stagnation  acoustic  speed,  and  r,c  is  the 
rotational  relaxation  time.  Figure  15  shows  the  freezing  parameter, 
d*  ■'  *tc  Be  ,  as  a  function  of  Mach  number  for  the  sonic  orifice  flow  from 
Ref.  10  and  for  two  nozzles  along  streamlines  bounding  0,  50  percent, 
and  90  percent  of  the  mass  flux.  The  rotational  freezing  parameter  can 
be  put  in  a  convenient  form  for  a  rigid  sphere  molecule: 


d*/r 


r 

e 


3.9  x  10* 


3.9  ,  .O'  (£-)  d- 


for  N2  :  -p—  is  the  supply  pressure  in  atm  and  d*  is  in  cm.  In  Figs.  2  and 
X  the  contours  of  rotational  freezing  points  are  superimposed  for  a  throat 
diameter  of  1.  0  mm.  For  the  nozzles  used  in  these  experiments,  rota¬ 
tional  freezing  effects  on  the  flow  field  may  be  anticipated  at  supply  pres¬ 
sures  less  than  about  2  atm. 


The  freezing  of  the  translational  mode  is  downstream  of  the  rotational 
freezing  point  and  will  not  influence  the  flow  field  in  terms  of  velocity  and 
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density.  Of  course  the  pressure  and  temperature  are  entirely  different 
when  translational  freezing  occurs  since  the  temperature  approaches  a 
constant  instead  of  zero.  This  phenomenon  has  been  widely  investi¬ 
gated  in  the  literature  of  molecular  beam  velocity  distributions  (for 
instance.  Ref.  11  is  a  good  example  ). 

2.1.5  Summary  of  Limitations  of  Applicability  of  Characteristics  Calculations 

1.  Condensation:  At  supply  pressures  from  1-7  atm  and  tempera¬ 
ture  of  300°K,  condensation  may  occur  for  nozzle  area 

ratio  >70  for  nitrogen.  For  argon,  there  are  insufficient  data 
to  predict  condensation,  cf  Appendix  I. 

2.  Viscosity:  The  nozzle  boundary  layer  reduces  the  exit  Mach  num¬ 
ber  and  the  density  in  the  outer  edge  of  the  jet.  The  effect  of 
viscosity  is  expected  to  be  negligible  outside  the  boundary  layer 
because  the  velocity  is  practically  uniform. 

3.  Continuum  Effects:  At  supply  conditions  where  the  product  of 
reservoir  density  and  throat  radius,  -jp-  r*  >  0.5  amagat-mm, 
the  flow  field  will  be  unchanged  by  the  transition  from  a  con¬ 
tinuous  medium, 

4.  Rotational  Freezing:  Rotational  freezing  will  not  affect  the  flow 
field  appreciably  if  r*  >  1  amagat-mm. 

SECTION  III 
APPARATUS 


3.1  GENERAL  TEST  DESCRIPTION 

The  tests  were  performed  in  the  Aerospace  Research  Chamber 
ARC  (8V).  The  nozzle  flow  was  pumped  with  20°K  cryosurfaces  and  a 
diffusion  pump-mechanical  pump  system.  The  jet  exhaust  flow  field  was 
investigated  with  flat  plates  nearly  normal  to  the  flow  direction  pitot 
v  probes,  and  streamline  vanes. 

3.2  AEROSPACE  RESEARCH  CHAMBER 

The  Aerospace  Research  Chamber  (8V)  is  a  nickel-plated  steel  cham¬ 
ber,  8  ft  in  diameter  by  16  ft  long:  a  removable  ellipsoidal  dished  head 
on  one  end  provides  access  to  the  chamber.  A  6-in.  diffusion  pump  in 
series  with  a  mechanical  pump  is  installed  on  the  chamber.  A  liquid- 
nitrogen-cooled  liner  shrouds  the  sides  and  downstream  end  and  provides 
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a  radiative  heat  sink  (Fig.  16).  A  1-kw  gaseous  helium  refrigerator  is 
available  for  cooling  240  ft^  of  20°K  cryosurfaces.  The  pumping  system 
is  able  to  maintain  10'4  to  1 0 ~ ^  torr  at  test  nozzle  mass  flow  rates  of 
up  to  1.0  gm/sec  with  argon  and  nitrogen.  Heference  12  gives  further 
information  on  the  chamber  design.  A  schematic  of  the  nozzle  installa¬ 
tion  is  shown  in  Fig.  16. 


3.3  NOZZLE  ASSEMBLY 

The  test  articles  were  a  15°  half-angle  conical  segmented  nozzle 
and  a  sonic  orifice  (Fig.  17).  The  throat  diameter  for  all  nozzles  and 
the  orifice  was  0.  044  inch.  Detachable  segments  provided  area  ratios 
of  3.  73:1,  12.8:1,  36.9:1,  69.6:1,  and  207:l. 

A  stainless  steel  stagnation  chamber  supported  the  nozzle  and  the 
propellant  supply  line.  The  stagnation  chamber  and  nozzle  assembly 
were  rotated  in  the  vertical  plane  through  50  degrees  of  pitch  from  the 
horizontal  position  by  a  10:1  worm  gear  driven  by  a  shaft  extending 
through  the  door.  This  shaft  also  provided  nozzle  assembly  axial  move¬ 
ment  of  20  inches  by  a  push-pull  motion  of  the  handle.  The  combination 
of  pitch  and  axial  movement  was  used  to  map  the  plumes  using  stationary 
instrumentation.  Nozzle  angle  was  measured  by  a  potentiometer 
attached  to  the  rotation  shaft;  axial  position  was  measured  externally  on 
the  driving  shaft  by  a  tape  measure. 


3.4  PROPELLANT  SUPPLY  SYSTEMS 

A  schematic  of  the  nozzle  propellant  supply  system  is  shown  in 
Fig.  18.  Stagnation  pressure  was  controlled  by  line  regulators  and 
monitored  by  a  Kollsman  gage. 

3.5  CANTILEVERED  FLAT  PLATE 

The  jet  properties  were  defined  by  measuring  the  forces  on 
1  in.  x  1  in.  and  1/4  in.  x  1/4  in.  square  flat  plates  nearly  normal  to  the 
flow  direction.  The  plates  were  maintained  stationary  while  the  nozzle 
was  positioned  as  described  in  Section  3.  3.  The  flat  plate  was  mounted 
on  a  36-in.  aluminum  beam.  The  lower  end  of  the  beam  was  clamped  to 
a  support  mount  in  the  chamber  floor.  The  beam  was  shielded  from  the 
flow  by  3- in.  diameter  tubing  and  aluminum  sheet  (Fig.  19).  Four 
120-ohm  thin-film  strain  gages  were  cemented  to  the  base  of  the  aluminum 
beam.  These  gages  yielded  an  output  proportional  to  the  horizontal  force 
on  the  flat  plate.  Eddy-current  damping  was  employed  to  diminish  the 
effect  of  vibration  on  strain- gage  output. 
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The  flat  plate  was  calibrated  after  installation  as  illustrated  in 
Fig.  20.  The  method  provided  a  horizontal  resultant  from  a  vertical 
force  by  the  use  of  a  45-deg  inclined  thread.  Horizontal  force  and 
plate  displacement  were  calibrated  with  strain  gage  output.  The  cali¬ 
bration  curve  is  shown  in  Fig.  21. 


3.6  PITOT  PROBES 

Pitot  pressure  was  determined  by  two  pitot  probes  {1/16  in.  OD, 

1/32  in.  ID)  attached  to  5  and  15-psi  full-scale  transducers,  respectively. 
A  swing  arm  supported  the  pitot  probe  transducer  units  and  allowed  the 
probes  to  be  swung  away  from  the  nozzle  during  flat  plate  and  stream¬ 
line  vane  runs.  Pitot  pressure  was  investigated  in  the  region  0  to  1-in. 
from  the  nozzle  exit.  Minimum  resolution  for  the  transducers  was 
approximately  0.  05  psi. 


3.7  STREAMLINE  VANE  ASSEMBLY 

The  streamline  vanes  were  constructed  of  0.  0015-in.  thick  steel 
shim  stock.  The  vanes  were  mounted  on  steel  pins  protruding  from  the 
top  and  bottom  of  the  assembly  which  permitted  the  vanes  to  rotate  with 
a  minimum  of  force  to  align  with  the  flow  direction  (Fig.  23).  Photo¬ 
graphs  of  the  vanes  were  made  through  a  port  in  the  top  of  the  chamber. 
Two  500-watt  flood  lamps  were  installed  above  the  vane  assembly  to 
illuminate  the  vanes. 

3.8  PRESSURE  AND  TEMPERATURE  INSTRUMENTATION 

Ambient  pressure  was  monitored  by  an  ionization  gage  mounted  on 
the  nozzle  assembly  and  pointed  toward  the  chamber  wall.  Copper- 
constantan  thermocouples,  referenced  to  liquid  nitrogen,  were  used  to 
monitor  temperatures  of  the  insulated  heater  box,  flat  plate  strain  gages, 
and  nozzle  surface.  A  thermocouple  probe  was  inserted  into  the  nozzle 
stagnation  chamber  to  monitor  propellant  total  temperature. 


3.9  DATA  RECORDING  SYSTEM 

Cantilevered  plate  forces,  nozzle  pitch  angles,  pressures,  and 
temperatures  were  recorded  on  a  recording  digital  voltmeter.  Stream¬ 
line  direction  vane  photographs  were  taken  with  a  7  0-mm  camera  using 
Tri-X  film. 
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SECTION  IV 
PROCEDURE 


4.1  FORCE  AND  PRESSURE  RUNS 

During  cooldown  of  the  AHC  (8V)  cryosurfaces  to  operating  tem¬ 
perature,  the  beam  temperature  near  the  strain  gages  was  monitored 
to  assure  that  no  zero  shift  occurred.  Propellant  flow  was  then 
started  and  continued  until  equilibrium  pressure  had  been  reached. 

The  nozzle  assembly  was  then  moved  to  the  most  remote  position  from 
the  flat  plate  and  was  rotated  upward  to  the  50- deg  position. 

The  nozzle  was  then  rotated  from  50  to  -50  deg  and  returned  while 
recording  the  force  on  the  plate.  This  was  done  at  a  series  of  different 
axial  positions.  For  the  1-in.  plate,  sweeps  were  made  at  2-in.  axial 
intervals  at  distances  between  8  and  29  inches  from  the  nozzle  exit  and 
for  the  1/4- in.  plate  at  1-in.  axial  intervals  at  distances  between  1  and 
10  inches  from  the  nozzle  exit.  Pitot  pressure  sweeps  were  made  at 
l/4-in.  axial  intervals  at  distances  0  to  1  inch  from  the  nozzle  exit, 
although  probe  contact  with  the  nozzle  prevented  full  ±50-degree  sweeps. 
Each  run  took  approximately  10  minutes. 


4.2  STREAMLINE  VANE  RUNS 

The  startup  procedures  were  similar  to  those  described  for  the  force 
runs.  When  flow  was  stabilized,  the  flood  lamps  were  energized  and  the 
nozzle  moved  to  the  26-in.  axial  position  and  set  at  zero  degrees. 
Photographs  were  then  made  of  the  vane  angular  positions  at  2 -in.  axial 
intervals  from  26  inches  to  6  inches. 


4.3  RANGE  OF  VARIABLES 

Nozzles  with  stagnation  pressures  of  35,  50  and  100  psia  and  a  stagna¬ 
tion  temperature  of  300°K  were  investigated.  .Argon  and  nitrogen  were 
used  as  nozzle  fluids.  Nozzle  area  ratios  of  1.  0.  3.  73,  12.  8,  36.  9,  69.  6 
and  207  were  tested. 


4.4  DATA  REDUCTION 


Nondimensionalized  flat  plate  force  and  position  data  were  reduced 
by  a  computer  program.  Plots  of  nondimensionalized  force  versus 
angular  position  for  various  radial  positions  from  the  nozzle  exit  were 
obtained  from  the  cantilevered  flat  plate  runs  (typical,  Fig.  23).  A 
fifth-order,  least-squares  curve  was  fitted  to  the  experimental  data, 


12 


A  E  DC-TR-65-151 


and  the  constant-force  contours  were  crossplotted  from  the  fitted  curves. 
Data  point  scatter  bands  were  determined  from  the  nondimensionalized 
force  versus  angular  position  plots  and  is  indicated  on  the  corresponding 
constant-force  maps. 


SECTION  V 
RESULTS 


Computed  Mach  number  distribution  and  boundary  layer  displace¬ 
ment  thickness  (Cohen  and  Reshotko  solution)  in  the  nozzles  for  y  =  1.  4 
and  y  =  1.  67  are  shown  (Figs.  14a  and  14b).  The  centerline  exit  plane 
Mach  number  obtained  from  pitot  probe  results  for  Ae/A*  =  207  are  com¬ 
pared  with  inviscid  and  viscous  Mach  number  solutions.  Experimental 
Mach  numbers  were  computed  for  isentropic  flow  using  the  Rayleigh 
formula  for  pitot  pressure. 

Experimental  constant  force  contours  and  centerline  force  plots  are 
compared  with  the  characteristics  method  calculations  for  Ae/A*=  1.  0 
through  207  with  nitrogen  at  Pc  =  100  psia  and  for  Ae/A*  =1.0  and  3.  73 
with  argon  at  Pc  =  100  psia  (Figs.  24  through  33).  The  force  plots  for 
Ae/A*  =  1.  0  with  argon  and  nitrogen  are  compared  with  the  predictions  of 
Sherman  (Ref.  13)  since  the  present  characteristics  program  could  not 
be  used  for  Me  =  1  (sonic  orifice)  at  the  exit.  Sherman  calculated  by  the 
method  of  characteristics  the  centerline  Mach  number  of  the  sonic  orifice 
exhaust  flow  field  by  assuming  an  exit  Mach  number  slightly  greater  than 
unity.  Density  and  Mach  number  values  were  obtained  from  the  charac¬ 
teristics  method  calculations;  normal  flat  plate  force  was  determined  by 
the  method  outlined  in  Appendix  II  and  the  solution  compared  with  ex¬ 
perimental  results.  Knudsen  number,  based  on  the  flat  plate  dimension 
(1/4  in.  and  1  in.  ),  is  shown  in  Fig.  33. 

The  random  error  in  data,  indicated  by  the  maximum  deviations  from 
the  least  squares  curve  fit,  is  shown  on  the  constant  force  data  points. 

The  systematic  error  is  summarized  in  Table  I.  The  tare  drag  on  the 
1/4  in.  x  1/4  in.  flat  plate  support  beyond  the  windshield  caused  the  1/4- in. 
flat  plate  force  readings  to  be  approximately  12  percent  high.  The  flat 
plate  calibration  (Fig.  21)  indicates  a  possible  zero  shift  or  nonlinearity 
of  the  balance  of  approximately  ±0.  5  x  10"^  lb.  The  nozzle  chamber  pres¬ 
sure  deviation  is  negligible  (±0.5  percent).  The  possible  systematic 
errors  irh  the  data  are  shown  on  the  constant  force  contours  and  centerline 
plots. 

Streamline  vane  photos  for  Ae/A*=  1. 0  through  207  with  nitrogen  and 
for  Ar/A*=  1.0  and  3.  73  with  argon  are  compared  with  characteristics 
predictions  (Figs.  37  through  44). 


13 


AEDC-TR -65-151 


SECTION  VI 

DISCUSSION  OF  RESULTS 


6.1  EXIT  MACH  NUMBER 

Theoretical  calculations  and  experimental  test  results  both  indicate 
that  boundary  layer  effects  on  the  exit  Mach  number  for  nitrogen  at 
Pc  =  100  psia  are  negligible.  The  actual  Mach  number  distribution  across 
the  exit  plane  could  not  be  determined  because  of  the  small  size  of  the 
nozzle.  The  trend  in  the  pitot  results  for  argon  is  inconclusive  because 
of  possible  condensation  in  all  nozzles.  Exit  Mach  numbers  for  argon 
could  not  be  determined  because  static  pressure  was  not  monitored  and 
cannot  be  inferred  from  Ppjt0t/Pc  since  isentropic  theory  is  not  appli¬ 
cable  if  condensation  shocks  occur  in  the  nozzle. 


6.2  EXHAUST  FLOW  FIELD  -  NITROGEN 

Constant  force  contours  for  n2  suggest  that  noncondensed  flow 
existed  for  all  area  ratio  nozzles  since  all  forces  fell  inside  the  theoreti¬ 
cal  maximum  expansion  angle.  The  flat  plate  results  cover  the  plume 
expansion  up  to  a  mean  free  path  of  approximately  5  cm  (Kn  =  2  for  1-in. 
flat  plate)  based  on  the  assumption  p  A  =  constant. 

The  1-in.  flat  plate  results  show  agreement  with  the  predictions 
except  near  the  axis  where  the  measured  density  is  lower  than  predicted. 
Computation  difficulties  with  the  method  of  characteristics  along  the 
centerline  may  account  for  part  of  this  deviation.  Characteristics  pre¬ 
dictions  in  this  region  were  approximated  by  fairing  from  the  outer 
portions  of  each  curve.  In  Fig.  29a,  the  (K/A)PC-  =  2.  5  x  10  ^  contour 
for  Pr.  =  35  psia  shows  the  trend  of  the  force  field  with  decreasing  nozzle 
chamber  pressure.  The  outward  shift  of  the  2.  5  x  10"^  contour  is  pre¬ 
dicted  by  the  outward  shift  of  the  theoretical  contour  with  increasing 
Knudsen  number.  The  experimental  contours  (  Pr  =  100,  35)  show  the 
same  trend  as  the  corresponding  theoretical  contours. 

The  1/4 -in,  flat  plate  results  show  higher  density  than  predicted 
near  the  axis  and  a  lower  density  than  predicted  in  the  outer  portions  of 
the  plume.  The  disagreement  in  the  outer  portion  could  have  been 
caused  by  the  boundary  layer  effects  and/or  by  the  force  gradient  on  the 
flat  plate.  The  flow  impingement  on  the  supporting  beam  would  result 
in  forces  that  are  approximately  12  percent  high  (see  Table  I). 

Axial  distribution  of  the  flat  plate  force  shows  the  effects  of  varying 
stagnation  pressure,  thereby  changing  the  Knudsen  number  and  drag 
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coefficient.  The  effect  is  in  the  expected  direction  with  the  lower  stagna¬ 
tion  pressure  yielding  a  higher  (F/A)PC  .  The  effect  of  plate  size  on  the 
Knudsen  number  (and  CD)  is  in  the  anticipated  direction,  although  the  dif¬ 
ference  in  the  results  cannot  entirely  be  accounted  for  by  the  effects  of 
plate  size  alone.  Even  if  tare  effects  are  taken  into  account,  the  1/4-in. 
flat  plate  results  (Ae/A*  =  3.  73 — ►  207)  indicate  a  higher  density  while  the 
1-in.  flat  plate  indicates  a  lower  density  than  predicted.  This  discrepancy 
could  be  due  to  a  systematic  error,  i.  e.  ,  a  zero  shift  or  a  nonlinearity 
which  is  larger  than  the  magnitude  assumed.  The  indicated  density  by  both 
flat  plates  for  the  sonic  orifice  is  higher  than  predicted  by  Sherman 
(Ref.  13). 

The  streamline  directions  indicated  by  the  vanes  align  with  the  pre¬ 
dicted  directions.  The  streamlines  become  linear,  as  in  source  flow, 
a  short  distance  from  the  nozzle  exit.  The  apparent  source  position 
varies  with  area  ratio.  The  streamline  vane  misalignment  in  the  outer 
portions  of  the  flow  field  is  large  because  of  low  momentum  flow  and  static 
moments  on  the  vanes. 

The  agreement  of  the  constant  force  contours  with  the  predictions  in 
the  outer,  near  free  molecule,  portion  of  the  plume  indicates  that  the 
method  of  characteristics  is  useful  for  free  expansions  from  initial  con¬ 
tinuum  conditions  into  a  rarefied  environment.  Flow  velocity  and  direc¬ 
tion  are  established  in  the  continuum  (near  the  exit)  flow  regime.  For 
y  ~  1*4,  99  percent  Vmax  is  reached  by  M  =  16  which  occurs  in  the  con¬ 
tinuum  region  of  the  plumes  investigated. 

Even  though  Mach  number  and  temperature  are  frozen  at  some 
point  in  the  field,  density  can  still  be  predicted  from  velocity  and  mass 
flow.  The  normal  flat  plate,  which  measures  momentum,  is  then  a  den¬ 
sity  indicator  since  velocity  is  essentially  constant.  The  similarity  of 
the  continuum  and  free  molecular  flow  relations  is  discussed  in  Sec¬ 
tion  2.1, 

6.3  EXHAUST  FLOW  FIELD -ARGON 

The  results  for  argon  are  presented,  even  though  no  direct  com¬ 
parison  with  characteristics  predictions  of  the  present  program  is  pos¬ 
sible.  The  sonic  orifice  axial  distribution  of  flat  plate  force  is  compared 
with  the  predictions  of  Sherman  (Ref.  13). 

The  constant  force  maps  for  Ac/A*  =  12.  8  to  207  show  expansion 
angles  greater  than  the  maximum  expansion  which  suggests  condensation; 
however,  Ae/A*  =1.0  and  3.  73  did  not.  Decreasing  the  stagnation  pres¬ 
sure  in  the  Ac/A*=  12.  8  nozzle  did  not  cause  a  significant  decrease  in  the 
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expansion;  this  indicates  that  the  condensation  shock  was  not  positioned 
near  the  exit.  Decreasing  the  stagnation  pressure  would  increase  the 
supersaturation  of  the  flow,  thereby  positioning  the  condensation  shock 
farther  downstream  in  the  nozzle.  The  similarity  of  the  Ae''A*  =  3.  73 
constant  force  maps  for  nitrogen  and  argon  should  be  noted.  The  simi¬ 
larity  was  not  anticipated  because  of  the  different  expansion  properties 
of  the  two  gases. 

The  axial  distribution  of  flat  plate  force  for  Ac./A'c  =1.0  indicates 
a  higher  density  for  both  plates  than  calculated  by  Sherman  (Ref.  13). 

A  slight  Knudsen  number  effect  is  indicated  with  the  changing  plate  size 
and  stagnation  pressure,  whereas  theoretical  calculations  indicated  no 
effect  on  the  Knudsen  number. 

The  streamlines  indicated  by  the  vanes  align  well  with  the  esti¬ 
mated  streamlines,  which  are  linear  extensions  of  the  calculated  stream¬ 
lines. 

Even  though  Mach  number  is  not  a  physically  significant  quantity  in 
free  expansion  analyses,  the  method  of  characteristics  relations  in¬ 
volving  Mach  number  still  yield  meaningful  density  predictions.  The 
flow  fields  tested  essentially  become  source-like  early  in  the  expansion; 

i.e.  ,  the  density  varies  inversely  with  the  square  of  the  distance. 


SECTION  V|| 
CONCLUSIONS 


The  analysis  of  a  flow  field  expanding  to  low  background  pressures 
and  the  experimental  results  have  demonstrated  that; 

1.  The  method  of  characteristics  yields  accurate  density  predictions 
for  expansions  from  continuum  to  the  rarefied  flow  regimes  if 

(a)  the  transition  to  free  molecular  flow  takes  place  in  the  source¬ 
like  flow  region  and  (b)  condensation  does  not  occur  in  the  con¬ 
tinuum  or  transition  region. 

2.  Condensation  in  the  nozzle  apparently  causes  expansion  greater 
than  predicted  for  argon. 

3.  The  normal  flat  plate  technique  is  a  useful  tool  for  rarefied  flow 
investigations  and  is  limited  primarily  by  force  balance  sensitivity. 

4.  The  light  streamline  vane  is  useful  for  determining  flow  direction 
in  the  continuum  and  transition  regimes. 
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TABLE  I 

SYSTEMATIC  ERROR  OF  FORCE  BALANCE 


Pc  (»*-.. 

f  AF(lb) 

AP{ 

A  l  (F/A) / Pc  ]  1/4"  Flat  Plate 

A[(F/A)/PC  3  1"  Flat  Plate 

100 

±.5  x  10-4 

±.57« 

±  .8  x  10“5 

+  .5  x  10-6 

50 

1.5  x  10“4 

+.5% 

±.16  x  10-5 

+1.0  x  10-6 

35 

+.5  x  10“4 

+.57) 

±2.4  x  10“5 

±1.5  x  10-6 

NOTE:  Additional  frontal  area  exposed  to  the  flow  for  the  1/4"  x  1/4”  flat  plate  by  the 
supporting  beam  is  0.0078  sq.  in.  Therefore,  the  1/4"  x  1/4"  flat  plate  readings 
are  estimated  to  be  12%  high.  (See  Fig.  17) 
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Fiq.  2  Method  of  Characteristics  Flow  Field  Solution;  A„/A*  =  3.73,  y  -  1.4,  M0  =  2.88 
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Fig.  5  Method  of  Characteristics  Flow  Field  Solution;  AB/A 
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Fig,  11  Method  of  Characteristics  Flow  Field  Solution;  Ae/A*  =  207,  y  =  1-67,  Me  14.9 
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Fig.  12  Effect  of  Nozxle  Exit  Area  Ratio;  Pt 
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Fig.  16  Nozzle  Installation  —ARC  (8V) 
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Fig.  18  Propellant  Supply  System 


Fig.  19  ConlilevtnJ  Flat  Plate 
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Fig.  20  Flot  Plate  Calibration 
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Fig.  21  Flat  Plat«  Calibration  Curve 
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Fig.  23  Force  Distribution  (Typical);  Nitrogen,  Pc  =  50  psia, 
A/A*  =  12.8,  y  =  1.4,  1- x  1-in.  Flat  Plate 


Fig.  24a  Constant  Force  Map;  A„/A*  =  1,  Nitrogen,  Pc  =  100  psia 
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Fig.  25a  Constant  Force  Map;  Ao/A*  =  3.73,  Nitrogen,  Pc  =  100  psio 
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Fig.  26a  Constant  Force  Map;  Ae/A*  =  12.1 
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Fig.  27a  Constant  Force 
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Fig.  28a  Axial  Distribution  of  Flat  Plate  Force;  Ae/A*  =  69.6,  Nitrogen,  P 


Fig.  28b  Axial  Distribution 
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Fig.  29a  Constant  Force  Map;  AB/A*  =  207,  Nitrogen 
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Fig.  31a  Constant  Force  Map,  Argon,  Ae/A*  1.0,  Pc  =  100  psia 
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31b  Axial  Distribution  of  Flat  Plato  F 
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Fig.  32a  Constant  Force  Map;  Ac/A*  =  3.7 
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Fig.  32b  Axial  Distribution  of  Flat  Plate  Force;  Ae/A*  =  3.73,  Argon 
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Fig.  34  Flow  Direction  by  Streamline  Vanes;  Ae/A*  =  1.0,  Nitrogen,  Pc  =  100  psio 
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Fig.  40  Flow  Direction  by  Streamline  Vanes;  Ae/A*  =  1.0,  Argon,  Pc 
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APPENDIX  I 

CONDENSATION  IN  SUPERSONIC  NOZZLE  FLOW 


When  an  ideal  gas  expands  is ent Topically  through  a  nozzle,  the  tem¬ 
perature  and  pressure  decrease  such  that  P/(T) y/y~  1  =  Const.  This 
isentrope  will  cross  the  saturated  vapor  pressure,  Pv  =  f(T)  ,  on  a 
P  -  T  plane  and  the  gas  will  tend  to  condense  in  liquid  or  solid  drop¬ 
lets.  The  amount  of  supersaturation  that  can  be  sustained  before  con¬ 
densation  is  observed  has  been  the  subject  of  a  great  deal  of  theoretical 
and  experimental  work.  The  experimental  data  from  various  wind  tun¬ 
nels  show  the  amount  of  supersaturation  that  can  be  tolerated  in  terms 
of  the  effects  on  certain  measurable  quantities  for  the  specific  nozzle. 
The  operating  limits  for  the  facilities  are  satisfactorily  defined.  How¬ 
ever,  the  generalization  of  these  results  to  predict  condensation  for 
arbitrary  conditions  is  not  possible.  The  classical  theory  of  condensa¬ 
tion  outlined  below  is  a  valid  description  of  the  phenomenon  at  higher 
pressures,  but  it  appears  to  give  ambiguous  results  in  the  usual  wind 
tunnel  operating  range. 


THE  CLASSICAL  CONDENSATION  THEORY 


References  1  and  2  discuss  the  theory  of  condensation  in  isentropic 
expansions  developed  by  Oswatitsch,  Volmer,  Becker  and  Doring, 
Frenkel  and  others  to  account  for  observations  in  supersonic  wind  tun¬ 
nels.  When  the  gas  is  in  a  supersaturated  state,  droplets  will  form  by 
the  aggregation  of  molecules  and  the  transfer  of  the  heat  condensation  by 
the  surrounding  gas.  The  size  and  rate  of  formation  of  these  droplets  is 
predicted  from  the  theory.  The  concept  of  the  critical  droplet  size  is 
essential  to  the  theory  since  it  is  the  formation  rate  and  growth  rate  of 
these  that  determine  condensation. 


The  critical  droplet  is  of  a  size  that  it  is  unstable  with  respect  to 
the  surrounding  saturated  vapor;  if  a  molecule  evaporates,  the  critical 
droplet  will  continue  to  shrink,  or  if  one  is  added  it  will  grow.  When 
droplets  of  critical  size  are  formed,  condensation  proceeds  rapidly. 
Both  the  critical  radius  and  rate  of  formation  of  droplets  of  critical  size 
depend  upon  the  state  of  the  gas  and  the  surface  tension  which  binds  a 
molecule  to  the  droplet: 


rcrit 


KT  In  p/p^ 


critical  radius 


l  4  P  .  1 2  n  cr  4  nr  a  ( r  CI:,  )* 

v  —  -  t  tt-  ir"li)  - — 

rate  of  condensation 
fraction 
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where:  c 

Vj 


T  &  P 


Surface  tension 

Volume  of  molecule  in  liquid 

Mass  of  molecule  in  liquid 

Equilibrium  vapor  pressure  at  T 

Temperature  and  pressure  of  supersaturated  vapor 


The  situation  on  a  P  -  T  plane  is  shown  schematically  in  the  figure 
below: 


Pressure 


Along  the  isentropic  expansion,  1  ,  the  critical  droplet  size  decreases, 
and  the  formation  rate  increases  after  the  saturation  line  is  crossed.  This 
rate  increases  slowly  at  first  and  then  rises  rapidly  as  indicated  by  the 

lines  —  ~  =  1  and  106  per  second.  Knowing  the  flow  velocity,  one  can 
calculate  the  fraction  of  the  flow  condensed  by  integrating  along  the  isen- 
trope.  Lowering  the  supply  pressure  or  increasing  stagnation  tempera¬ 
ture  shifts  the  isentrope  from  I,  to  Ia  and  evidently  increases  the  amount 
of  super  saturation.  It  is  also  evident  that  the  critical  droplet  radius  is  de¬ 
creased  at  the  point  where  appreciable  formation  rates  are  encountered. 

The  last  observation  and  the  absence  of  any  evidence  of  condensation 
in  the  NOL  wind  tunnel  lead  Bogdonoff  and  Lees  (Ref.  1)  to  postulate  a 
breakdown  of  the  theory  when  the  droplet  radius  is  small,  i.  e.  ,  of  the 
order  of  the  gas  molecule  radius.  For  small  droplets  the  surface  tension 
is  reduced  from  that  in  an  infinite  body  of  liquid  because  fewer  molecules 
are  near  a  surface  molecule  on  a  surface  of  small  radius.  This  leads  to 
the  conclusion  that  droplets  predicted  by  the  theory  are  smaller  than 
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critical  size  and  therefore  they  evaporate  rapidly  instead  of  growing. 

A  cut-off  in  the  nucleation  rate  will  appear  when  the  predicted  droplet 
radius  is  small  as  indicated  on  the  sketch  above.  Therefore,  when  an 
isentrope  intersects  the  saturation  curve  below  a  fixed  temperature, 
condensation  will  not  occur  at  all. 

A  modification  of  the  classical  theory  is  given  by  Stever  and 
Rathbun  (Ref.  2)  in  which  the  effect  of  droplet  radius  on  surface  tension 
is  computed.  The  result  of  this  refinement  is  to  shift  both  the  lines  of 
constant  critical  droplet  radius  and  condensation  rate  toward  the  satura¬ 
tion  line.  Therefore,  the  amount  of  supersaturation  estimated  is  con¬ 
siderably  less  than  the  original  theory  for  constant  surface  tension  would 
estimate.  Experimental  results  were  presented  which  tended  to  sub¬ 
stantiate  the  modified  theory,  and  it  was  suggested  that  the  observations 
of  Ref.  1  were  insensitive  to  the  existence  of  condensation.  However, 
the  breakdown  of  the  theory  for  small  droplets  is  not  precluded,  and  two 
of  the  data  points  appear  to  substantiate  it. 


EXPERIMENTAL  CONDENSATION  RESULTS 

The  most  recent  compilation  of  wind  tunnel  data  concerning  con¬ 
densation  are  from  Daum  (Ref.  3).  Figure  1-1  shows  the  data  of  this 
report  and  some  additional  points  from  AEDC  wind  tunnel  experience. 
Isentropic  lines,  the  condensation  rate,  and  critical  droplet  contours 
from  the  theory  of  Ref.  2  are  also  included.  These  data  support  the 
nucleation  theory  of  Ref.  2  as  well  as  the  hypothesis  of  Bogdonoff  and 
Lees,  at  least  to  the  extent  of  indicating  much  higher  supersaturations 
when  the  isentrope  crosses  the  saturation  curve  lower  than  10  mm  Hg 
pressure.  This  occurs  when  the  droplet  radius  is  less  than  about 
3  nitrogen  molecule  radii  at  high  nucleation  rates.  It  is  interesting  to 
observe  that  the  limit  of  the  data  falls  along  a  single  isentrope  which  is 
a  little  below  the  point  where  the  theory  appears  to  be  invalid.  The  lack 
of  data  below  this  limiting  isentrope  is  the  result  of  (1)  limitation  of  wind 
tunnel  operating  characteristics  or  (2)  densities  too  low  to  make  reliable 
measurements. 

The  quantitative  correlation  of  wind  tunnel  data  is  rendered  difficult 
by  the  variety  of  methods  used  to  detect  the  onset  of  condensation.  In 
Ref.  1  a  condensation  shock  (or  absence  of  it)  was  anticipated  to  be  ob¬ 
servable  by  the  schlieren  system.  In  Ref.  2  condensation  effects  were 
observed  (1)  in  the  nozzle  wall  static  pressure,  (2)  in  the  shock  wave 
angle  on  a  5.5°  wedge,  and  (3)  in  the  dispersion  of  light  by  droplets  in  the 
flow.  The  deviation  of  static  and  pitot  pressures  as  the  total  temperature 
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was  lowered  was  used  in  Ref.  3  to  indicate  the  onset  of  condensation. 
Reference  4  used  the  nozzle  wall  static  pressure  and  observed  no  meas¬ 
urable  effect  on  pitot  pressure.  In  Ref.  5  the  deviation  of  the  surface 
pressure  on  an  ogive- cylinder  was  observed.  Reference  6  employs  the 
test  section  static  pressure  derived  from  an  extrapolation  of  flat  plate 
induced  pressures  to  an  infinite  distance  downstream.  All  the  methods 
are  indirect  measurements  and  are  therefore  more  or  less  related  to 
condensation  itself.  Further,  the  regime  of  interest  in  the  present  ex¬ 
periment  (Tc:  =  300°K,  Pc  =  1  -  7  atm )  lies  in  the  area  of  condensation 
where  the  orthodox  nucleation  theory  becomes  doubtful. 


A  KINETIC  MODEL  METHOD  OF  CORRELATING  CONDENSATION  DATA 


Since  there  is  no  reliable  guide  in  the  conventional  nucleation  theory, 
a  simple  kinetic  model  appears  to  be  justified;  this  is  to  count  the  num¬ 
ber  of  molecular  collisions  in  an  isentropic  expansion  from  the  point  of 
saturation  to  condensation  in  order  to  correlate  the  experimental  data. 
The  collision  frequency  can  be  calculated  using  a  rigid  sphere  model: 


(IN 

dl 

V 

A. 

P 

T 

Sub  o 
Sub  c 


v  _  V.  r  _Ls  (P_\  I 

A  A0  po  yj  T 0  UJV 


JL 

Tr 


Mean  molecular  velocity 
Mean  free  path 

Density 
Temperature 
Standard  T  and  P 
Nozzle  reservoir 


The  collisions  per  unit  distance  traveled: 


rfN 

dx 


1  dN 
U  dt 


I  y-i 
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U  = 


Stream  speed 

Reservoir  acoustic  speed 
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ar.  =  Limit  speed 
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Then, 


is  the  number  of  collisions  between  and  x3  .  K 


is  a  constant  for  the  gas.  Since  the  quantity  inside  the  integral  is  a 
function  of  the  nozzle  area  ratio,  it  is  convenient  to  transform  the  equa¬ 
tion  by  taking  an  area  ratio-distance  relation  of  the  following  form: 


A_ 

A* 


+  const 


and  defining  the  constant  such  that  jC,„  corresponds  to  an  area  ratio  of  10. 
Then  the  number  of  collisions  between  an  area  ratio  of  10  and  is 

A 


/ 


f  /  A  \  <i  (A/A*  ) 

u^-r" 


f  {  — _ \  —  tl  m  a  x  P  ^ 

where  \  a*  /  u  pc  tc 


The  exponent,  n  is  nearly  unity  for  two-dimensional  nozzles,  two 
for  conical,  and  generally  larger  for  the  free  expansion  of  a  jet  into  a 
vacuum.  Figure  1-2  shows  the  effects  of  y  and  n  on  the  collision  num¬ 
ber  for  several  examples. 

The  theory  above  is  at  least  qualitatively  in  agreement  with  observa¬ 
tion.  Equation (4)  shows  that  a  reduction  of  the  reservoir  density  -p2- 
results  in  greater  supersaturation.  The  linear  effect  of  nozzle  size  (410) 
is  evident  in  comparing  the  greater  degree  of  supersaturation  (~  25°K) 
observed  in  a  small  nozzle  (-13  cm  long)  in  Ref.  4  compared  to  Ref.  5 
(-  6  m  long  nozzle)  where  it  was  only  a  degree  or  two.  It  is  interesting 
to  observe  that  the  integral  has  a  finite  limit;  therefore,  expansions  from 
a  low  density  reservoir  or  a  small  nozzle  could  be  expected  to  continue 
to  infinity  without  the  occurrence  of  condensation.  As  might  be  antici¬ 
pated,  the  more  rapid  the  expansion,  (n  >  1)  ,  the  fewer  are  the  col¬ 
lisions  and  therefore  greater  supersaturation  should  occur.  The  reason 
for  a  smaller  number  of  collisions  with  larger  ratio  of  specific  heats,  y , 

is  that  for  a  fixed  expansion  (jfc]  the  temperature  drops  further. 

To  correlate  condensation  results,  it  is  necessary  to  compute  the 
area  ratios  at  saturation  and  at  condensation  and  take  the  difference  in 
AN. 
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ESTIMATE  OF  THE  PROBABILITY  OF  CONDENSATION  IN  THE  PRESENT  EXPERIMENTS 

In  the  present  experiments,  conical  nozzles  with  area  ratios  up  to 
200  were  used.  Argon  and  nitrogen  were  the  test  gases  at  a  reservoir 
temperature  about  300°K  and  pressures  from  0.  7  to  7  atm.  It  was  ob¬ 
served  in  a  previous  test  (Ref.  AEDC-TDR-64-124)  that  the  exhaust 
plume  of  nitrogen  was  approximately  the  shape  calculated  by  the  method 
of  characteristics  whereas  argon  had  a  much  wider  expansion.  It  was 


suggested  that  a  condensation  shock  reduced  the  exit  Mach  number  in 

argon,  but  it  was 

absent  in  nitrogen.  The  classical 

nucleation  theory 

gives  the  following  results  for 

A/A*  =  70: 

Argon 

Nitrogen 

rc  ril 

1. 09  x  10"8 

cm  2.  34  x  10“8  cm 

Critical  Radius 

l_  dN 

N  dt 

2.  26  x  104/ 

sec  8.  95  x  104/sec 

Rate  of  Condensation 

Ax 

0.  13  cm 

0.  67  cm 

Distance  from  Satu¬ 

'  <_ 

ration  to  Exit 

Tc/M  1 

2.  74 

3.  27 

Proportional  to 

Tranart  Time 

l/r/  jUv *, 

1  /dN  \  Ax 

1.  44  x  104 

3. 28  x  104 

Proportional  to 
Condensation 

N  '  d  t  /  V'  Tc/M 

f>f  —  Mot.  us/. 

These  figures  indicate  practically  the  same  degree  of  condensation 
for  both  gases  and  give  no  clue  to  the  divergent  results  observed.  How¬ 
ever,  the  critical  radius  is  on  the  same  order  as  that  of  the  molecule  and 
this  indicates  the  nucleation  theory  is  not  applicable  in  the  regime  of  these 
experiments. 

Recourse  was  then  made  to  the  collision  theory  of  the  previous  sec¬ 
tion  by  correlating  the  data  of  Daum  (Ref.  3)  for  nitrogen.  Figure  1-3 
shows  the  number  of  collisions  for  nitrogen  and  argon  in  nozzles  with 
area  ratios  up  to  200.  The  data  points  of  Daum  and  Ref.  6  are  tabulated 
below  in  terms  of  number  of  collisions  between  saturation  and  condensation: 
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If  to  be  conservative  one  disregards  the  last  two  points,  the  super¬ 
saturation  temperature  increment  ATS  -  e  can  be  plotted  as  a  function  of 
the  number  of  collisions  ANS  _  e  ,  which  gives  the  following  estimates  for 
the  test  nozzles  with  nitrogen: 


A/A"  =  70  A-'A"  =  200 


An 

AT 

AT  |j  Air 

l  AT 

i  AT 

<°K) 

K 

(Calf) 

i  j  FjX  pariesjon )  |  K 

fCalO 

I.IjJC  pan*  ion 
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>  30’K 
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22 

22  1  ,  2. 0  x  10'2 

16 

34 

300 

2.  65  x  10~2 

1  15 

24  3.  36  x  in'2 

12 

l  36 

For  the  test  conditions  from  2  to  7  atm  supply  pressure,  condensa¬ 
tion  is  not  expected  for  the  70:1  area  ratio  nozzle  but  may  occur  at  the 
higher  pressure  in  the  200:1  nozzle 

The  case  for  argon  is  inconclusive  because  of  a  lack  of  experimental 
data  of  the  type  for  air  and  nitrogen.  A  comparison  between  nitrogen  and 
argon  is  made  below. 

AN/Ki,, 


Pc 

A/A*  -  70 

A/A*  =  200 

N, 

1 

0. 25  x  10-2 

0.  44  x  10"2 

10 

4.  4  x  1CT2 

6.  4  x  10“  2 

A 

1 

1.  73  x  1CT2 

1. 88  x  10"2 

10 

20.  9  x  IQ'2 

22. 4  x  IQ'2 

Although  argon  undergoes  4  to  6  times  as  many  collisions  at  the  same 
density,  the  decrease  in  density  of  argon  by  a  factor  of  10  did  not  alter  the 
shape  of  the  jet  expansion.  The  number  of  collisions  for  condensation  in 
a  monatomic  gas  is  evidently  different  from  that  for  a  diatomic  one.  Sys¬ 
tematic  experimentation  is  required  to  provide  consistent  results  for 
correlation  by  this  method. 
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Fig.  I-]  Summary  of  Condensation  Data 
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APPENDIX  II 

DRAG  OF  NORMAL  FLAT  PLATE  FROM  CONTINUUM  TO 
FREE-MOLECULE  FLOW  CONDITIONS 


Correlation  of  flow  characteristics  in  the  jet  exhaust  through  meas¬ 
urement  of  the  force  on  a  plate  normal  to  the  flow  requires  an  estimate 
of  the  drag- coefficient  over  a  wide  range  of  flow  regimes.  Actual  ex¬ 
perimental  data  in  the  transition  between  the  continuum  model  and  the 
free-molecule  regime  are  very  limited,  and  consistent  theoretical 
methods  are  lacking.  The  method  of  estimating  the  drag  coefficient  is  to 
assemble  experimental  data  giving  the  trends  with  departure  from  well 
established  limits  of  the  continuum  or  free  molecule  flow  regimes. 


CONTINUUM  LIMIT 

Hoerner  (Ref.  1)  summarizes  the  high  Reynolds  number  regime  drag 
data  obtained  from  wind  tunnels  and  ballistic  ranges.  Although  these 
data  were  obtained  from  discs  and  short  cylinders,  they  should  be  appli¬ 
cable  in  the  hypersonic  regime  of  interest  because  both  skin  friction  and 
base  drag  components  are  negligible.  In  Fig.  II- 1  the  base  drag  (M  <  4  ) 
is  subtracted  out  and  plotted  separately. 


FREE  MOLECULE  LIMIT 


The  free-molecule  formulas  of  Ref.  2  can  be  written  for  the  flat 
plate  with  diffuse  reflection  at  the  wall  temperature,  Tw  as: 


CD  =  ^ 

V  77  s 


_L  iJLZ. 
2  V 


Y  -  1  V  Tw/'I’t 

V  “’fell 


("  l  *  ‘  +  a/JL  [Hi  V  Ltl/It  -1  ( 

L  **  2  \  y  (U/UmBIj  J 


1  -  erfs  ) 


where 


u 


V  2  R  T 


=  1 1-  \]  =  Speed  ratio 
\  2 
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The  substitution, 

T  w  _  y  ~  1  s2  Tk 

I  x  y  (  II  /  U  m  a  it  )  I  T 

has  been  made  for  convenience  in  computation.  The  adiabatic  wall  condi¬ 
tion,  Taw  =  r  Ty  ,  is  shown  in  Fig.  II-I  for  the  recovery  factor. 


y  =  l  .4  and  f(s)  is  given  in  Ref.  2. 


TRANSITION  REGIME  (LOW  REYNOLDS  NO.  VISCOUS  EFFECTS) 

In  the  supersonic  low  density  regime,  the  effects  of  viscosity  (low 
Reynolds  number)  dominate  the  flow  field  about  a  blunt  body;  and  it  is 
assumed  that  pitot  pressure  gives  a  good  measure  of  the  variation  away 
from  the  continuum  drag  value.  Therefore,  it  is  assumed  that  the  ratio 
of  drag  coefficient  to  the  high  Reynolds  number  Cr  is  the  same  as  the 
ratio  of  pitot  pressure  to  the  Rayleigh  stagnation  pressure.  The  data  of 
Refs.  3,  4,  and  5  are  cross  plotted  to  provide  the  Mach  number  varia¬ 
tion  at  constant  Reynolds  numbers  from  1.  0  to  10.  The  Knudsen  number 

|^Kn  =  is  uniquely  related  to  Mach  and  Reynolds  numbers  by 

Kn  _  for  a  rigid  sphere  molecule  model.  Constant  Knudsen 

number  contours  are  included  in  Fig.  II- 1  to  provide  continuity  into  the 
hypersonic  regime. 

TRANSITION  REGIME  (HYPERSONIC  NEAR-FREE  MOLECULE  FLOW) 

Reference  6  provides  experimental  two-dimensional  flat  plate  drag 
data  at  Mach  numbers  from  6  to  10  and  Knudsen  numbers  from  about 
unity  to  20  or  30,  where  a  close  approach  to  the  free  molecule  value  is 
demonstrated.  Over  the  Mach  number  range  of  the  experiments,  the 
drag  coefficient  depended  only  on  Knudsen  number.  All  the  experimental 
data  are  summarized  for  the  insulated  wall  case  in  Fig.  II- 1  as  solid  lines. 
This  conforms  to  the  conditions  of  the  experimental  procedure  in  which 
the  supply  temperature  and  surrounding  chamber  wall  were  room 
temperature. 

DRAG  AT  EQUILIBRIUM  WALL  TEMPERATURE 

The  flat  plate  in  the  Aerospace  Research  Chamber  (8V)  is  exposed, 
in  addition  to  the  convective  heat  transfer  from  the  gas  stream,  to 
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radiation  from  the  liquid-nitrogen-cooled  shroud  and  cryopump  down¬ 
stream  and  to  the  room  temperature  chamber  end  upstream.  The 
radiative  heat  transfer  is 


qr„<!  =  HV  ~  TV)  ►  ^  (Tn4  -  l>‘) 

The  convective  rate  is 

iconv  ■=  s>  PV(-P  (Taw  -  T„) 

Equating  the  two  yields  a  fourth- order  equation: 


T 


4 


1’ 


W 


rF*  <  T»'  .  /«. 

- 5 - j'1'' 


If  the  speed  ratio  is  large,  the  accommodation  coefficients  and 
emissivity  (<)  are  unity,  and  y  -  1.4: 

s,  =  1/2  ^  y  1  =  0.  86  -  Stanton  Number 

k  =  5,67  x  10"  12  w/cm^/ (°K)^  -  Boltzman  Constant 


Cp  =  .24  x  4.  186  =  1.  005  joule/gm/°K  -  specific  heat 
V  =  7.40  x  104  v'Tt7300  -  fiow  velocity 

p  s  1.  14  x  10“^  Pt,  )  (p/pc)  gm/cc  -  density 


r 


T  F 


Tr 


S,  pv Cp 

21t 


1.  166  -  recovery  factor 
300°K  -  front  wall  temperature 
80°K  -  rear  wall  temperature 

6.60  k  10“  Pc  (p/p,.)  (CK)’ 


Pc  -  Supply  pressure  (atm) 

The  two  limits,  (a)  negligible  convective  heat  transfer  -  low  density  and 
(b)  negligible  radiative  heat  transfer  -  high  density,  give  wall  tempera¬ 
tures  of  2  50DK  and  350°K,  respectively.  The  limits  of  hypersonic  drag 
coefficient  are  2.  85  and  3.  0  for  these  cases.  It  is  evident  that  the 
equilibrium  wall  Temperature  depends  upon  the  product  Pc  (p/Pc)  .  and 
Fig.  II- 2  gives  the  variation  of  drag  coefficient  for  a  supply  temperature 
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of  300°K.  This  figure  can  be  used  to  correct  the  insulated  wall  flat 
plate  drag  for  the  conditions  of  the  test  at  large  speed  ratios.  It  is  prob¬ 
able  that  the  same  trend  will  be  exhibited  also  at  lower  speed  ratios  (or 
Mach  numbers).  In  Fig.  II- 3  the  insulated  flat  plate  drag  of  Fig.  II- 1 
are  faired  into  smooth  curves. 


REFERENCES 


1.  Hoerner,  S.  F.  Fluid  Dynamic  Drag;  Practical  Information  on 

Aerodynamic  Drag  and  Hydrodynamic  Resistance. 

Midland  Park,  N.  J.,  1958.  (Second  Edition) 

2.  Schaff,  S.  A.  and  Talbot,  L.  "Handbook  of  Supersonic  Aerodynamics- 

Mechanics  of  Rarefied  Gases . 11  NAVORD  Report  No.  1488, 
Vol.  5,  Section  16,  February  1959. 

3.  Enkenhus,  K.  R.  "Pressure  Probes  at  Very  Low  Density."  UTIA 

Report  No.  R-43  (AD  126534),  January  1957. 

4.  Kosterin,  S.  I.,  et  al.  "An  Investigation  into  the  Effect  of  the 

Rarefaction  of  a  Supersonic  Stream  on  the  Readings  of  a  Total- 
Head  Probe."  Translation,  International  Chem.  Engr. ,  Vol.  3, 
No.  2,  April  1963,  pp,  251-255. 

5.  Potter,  J.  L.  and  Bailey,  A.  B.  "Pressures  in  the  Stagnation 

Regions  of  Blunt  Bodies  in  the  Viscous-Layer  to  Merged- 
Layer  Regimes  of  Rarefied  Flow. "  AEDC-TDR-63-168,  1963. 

6.  Ko,  D.  R.  "Drag  of  a  Two-Dimensional  Strip  Normal  to  Hypersonic 

Near  Free  Molecule  Flow."  University  of  California  Report 
No.  AS-64-4,  January  1964. 


85 


Drag  Coefficient 


Ref.  1  -  Flat  Plate  and  Short  Cyl. 

Ref.  3  -  Pitot  Pressure 

Ref.  4  -  Pitot  Pressure 

Ref.  5  -  Pitot  Pressure 

Ref.  6  -  2-dimensional  Flat  Plate 
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Fig.  11-1  Flat  Plata  Drag  Coefficient  Experimental  Data 
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